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Abstract. The impedance measurement technique consists in that the phase-dependent (parametric) in- 
ductance of the system is probed by the classical tank circuit via measuring the voltage. The notion of the 
parametric inductance for the impedance measurement technique is revisited for the case when a quantum 
system is probed. Measurement of the quantum state of the system of superconducting circuits (qubits) is 
studied theoretically. It is shown that the result of the measurement is defined by the partial energy levels 
populations in the qubits. 
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1 Introduction 

The supercurrent /, flowing through a weak link between 
two bulk superconductors with the phase difference (j>, has 
the properties as a nonlinear inductor. This can be de- 
scribed by introducing the phase-dependent (parametric) 
inductance £ = (^o/27r) (dl/dcf>)~ . If the weak link is 
included in a ring, then the phase <p is related to the mag- 
netic flux <P, piercing the ring, <j> — 2tt<P / (Pq . The paramet- 
ric inductance can be measured [Tj and being inductively 
coupled to the resonant LC tank circuit provides the tool 
to measure the flux •P [2], [3]. The effective inductance 
of the tank circuit depends on the parametric inductance 
C(4>) and current I(<p) and thus the measurement in the 
tank circuit can be used for finding the inductance C [4], 
which is the so-called impedance measurement technique. 

The impedance measurement technique was recently 
applied for the measurement of the small currents in meso- 
scopic samples [5] and it was proposed for the descrip- 
tion of the currents in superconducting qubits [6], [7], [8], 
[9] and the series of the experimental results were ob- 
tained [TO] . However the theoretical works in this field 
consider mostly the ground state. If the superconduct- 
ing qubit is excited to the upper state, then the current 
in it has the probabilistic character, and in this way the 
parametric inductance depends not only on the clockwise 
(counter-clockwise) current value but also on the prob- 
abilities of the respective states. This consideration was 
used for the description of the phase-biased charge qubit 
pTj . And in this paper we study in detail the specifics 
of the impedance measurement technique when the quan- 
tum system is probed. The detailed presentation is aimed 
to show how the tank circuit is influenced by the para- 
metric inductance, and how this inductance have to be 



treated for the system of coupled superconducting circuits 
(qubits). For concreteness we consider the superconduct- 
ing circuits to be either flux or phase- biased charge qubits. 
The flux qubit [TJ] consists of a loop with three Josephson 
junctions. The phase-biased charge qubit [6], [7] consists 
of a loop with two closely situated Josephson junctions 
and with the gate which controls the charge on the island 
between the junctions. 



2 Tank circuit coupled to quantum object 
2.1 Equations for tank circuit 

The quantum system (coupled superconducting qubits) is 
considered to be weakly coupled via a mutual inductance 
M to the classical tank circuit. The circuit consists of the 
inductor Lt, capacitor Ct, and the resistor Rt connected 
in parallel. The tank circuit is biased by the current Iuas , 
and the voltage on it Vt can be measured. To obtain the 
equation for the voltage, we write down the system of 
equations, for the current in the three branches, namely, 
through the inductor (II), the capacitor (Ic), and the 
resistor (In) (see e.g. in the Chapter 14 of Ref. [5]): 



hias — II + Ic + Ir, 

Ic = e, e = CtVt, 
In, = Vt/Rt, 
V t = L T I L - $ e , 



(1) 
(2) 
(3) 
(4) 



Send offprint requests to: 



where e is the charge at the capacitor plate, the dot stands 
for the time derivative, is the flux through the tank 
circuit. This flux is the response of the quantum system 
to the flux, induced in it by the current II, and its time 
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derivative equals (see below for details): 

$e=Li L , (5) 

and thus equation (0| can be rewritten by introducing the 
effective inductance of the tank circuit L e f / : 



V T - L eff I L , 
L eff =Lt-L. 



(6) 
(7) 



Then from the system of equations (JTJ(4j) we derive the 
equation for the voltage in the tank circuit: 



C T V T + R^Vt + L~} f V T = I, 



'-bias • 



(8) 



2.2 Effective inductance of qubits 

Now we derive the relation ([5]); consider the flux <L> e = 
J2i , where is the flux induced by z-th qubit in 
the tank circuit: — M^I^ . Here MiT is the mutual 

(i) 

inductance of the qubit and the circuit, I qb is the current 
in the i-th qubit which equals to the expectation value 
of the current operator: lf£ = = Sp(j)I^j, where 

p is the reduced density matrix of the system of qubits. 
(Note that for one qubit, substituting <P e — MI q b, Eq. ((U 
coincides with Eq. (15) in Ref. [8].) 

The total flux that threads the loop of the i-th qubit 
#M consists of the external magnetic flux and the self- 
induced flux —Lilty (Li is the geometrical inductance of 
the loop): 

$M=*®-L i I$>. (9) 

This equation can be rewritten by introducing the para- 
metric inductance, 



r -l = 01 gb 



(10) 



to relate the variations of the external flux through the 
qubit 5<[>x and of the current in it 81 qb , as following: 

W W = s$® + = (A + Li) 51$. (11) 

Thus, we obtain the variation of the flux induced by the 
qubit in the tank circuit: 



The flux <P X in the i-th qubit is considered to consist 
of the fluxes induced by the tank circuit, M^Il, by the 
microwave source, <?ac smut, and by additional lines and 
by other qubits, 



<5<2>« 



(12) 



ifi 



(t) = M lT I L (t) + ft + #W sin cot. (13) 



Now let us recall that we consider the variation of the 
flux in order to calculate the derivative in time <? e which 
have to be substituted in Eq. ((4]) , which describe the tank 
circuit. We note that usually the dynamics of the tank 
circuit (with the frequency u> r f close to the resonant fre- 

quency lot = (LtCt) ) is significantly slower than the 
dynamics of a qubit driven by the microwave source at fre- 
quency to 3> lot- Hence, being interested in the response of 
the measurement system (that is of the tank circuit) , we 
average equations over the period 2tt/oo. After this averag- 

(i) 

ing the component <P a c smoot tends to zero and we have: 
5<I>x ~ Mix5lL- Here and below the time-averaging is as- 
sumed. (Note that here it is assumed that the expectation 

value for the current I^ b in Eq. (| 12(1 weakly depend on 
time during the time interval of the order of 2ix juo and its 
time dependence is defined by the tank circuit dynamics 
only.) Then it follows 



5<P 



St 



E, 



M?. 



(14) 



and Eq. ([5]) is obtained with the inductance L, which de- 
scribes the response of the quantum system to the tank 
circuit signal, given by: 



1 = 



Ci 



(15) 



Thus, we have obtained the system of equations (|7|), 
(0, (fT0| . and (|15p. which describe the interaction of the 
classical tank circuit and the quantum circuits (qubits). 
More accurate (quantum-mechanical) analysis would start 
the description from the Hamiltonian of the whole system 
in terms of the operators for both the tank circuit and 
the qubits with averaging the equations afterwards, as in 
Ref. [14] (see also the discussion about similar systems in 
[15] and [E]). Since this analysis would yield the same 

equations for the observable values (Vr — (Vr^, etc.), we 

do not consider this procedure here in detail. 



3 Analysis of the response of the tank circuit 

The measurement consists in biasing the tank circuit with 
the current Iuas = 1 A cos io r ft and measuring both the 
phase shift a and amplitude Va of the voltage Vr = 
Va cos(oo r ft + a). (The oscillations can be considered close 
to the harmonic form due to the small losses in the high- 
quality tank circuit which is weakly coupled to the non- 
linear qubits' inductances, see ([19]) and 1|20|) below.) Sub- 
stituting these expressions for Ibi as and Vr in Eq. {SJ and 
equating coefficients before sin ui r ft and cos co r ft, we ob- 
tain: 



tana 
V A 



Rt 



J eff 



rf j- -1 



LU r f \ LJ 

RtIa cos a. 



(16) 
(17) 
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Expression for the phase shift is simplified for the tank qubit as following (we rewrite the equations derived above): 



circuit driven at resonance, u) r f = lot = (ItC't) 



-1/2. 



tan a = Q 



L 



Q- 



L 



Lt - L ' L T 

E. k\L i L i ■ 
7. 



(18) 



Here it was assumed that L <C Lt and Li <ti £f, the latter 
inequality assumes that qubits' loops have small induc- 
tances, the former inequality is justified for large quality 
factor and small tank circuit-qubit coupling constants: 



Q = utRtCt > 1, 



Mir 
L^Lt 



«1, 



(19) 
(20) 



and also k 2 Q <1 is assumed. 

Alternatively to measuring the phase shift a, Eq. (fT8|) . 
the qubits' parametric inductances can be probed by mea- 
suring the amplitude Va of the voltage, which is optimal 
at the frequency ui defined by the relation (wy — ui) /lot = 
(2Q)" 1 , that is at ui = lj t (1 - (2Q)" 1 ), then from Eqs. 
1161)- (Eh it follows: 




-1/2 



(21) 



This relation can be rewritten, taking into account Eq. 
(|18p and assuming tan a < 1 (which is usually the case in 
experiment [TO]) in the form: 



Va 



-j=RtIa 



1 

1 + — tana| ( 



) r f —LOT 



(22) 



which shows the equivalence of the measurements via the 
amplitude and the phase shift of the tank circuit voltage. 
Since in practice it is more convenient to probe flux qubits 
via the phase shift [TO] , [T7] , we will consider in what fol- 
lows the phase shift only. 

For the case of small geometrical inductances Li, we 
can neglect the shielding current, then « 4^; we also 
assume the weak coupling of the qubits to the tank circuit 
(Eq. (|2T)|) ) and neglect the first term in Eq. fT3"]) , and define 
^dl as * ne corLS tant part of the flux through i-th qubit 
(in practice it changes adiabatically slow); hence 5<L>^ « 



(i) 



6<P d 1 } and for calculations Eq. ([T8]) is supplemented 



by the relation: 



U1 qb 



(23) 



Note that this value (namely, the r.h.s. of Eq. (|23|)) can 
also be interpreted as the magnetic susceptibility [pj , [T5] . 



4 Inductance of superconducting qubits 

Consider the case of a single superconducting qubit in 
more detail. The phase shift a probes the current in the 



tana 

Iqb 



k 2 QL£~ 

Olqb 

I) 



Sp (pi) 



(24) 
(25) 

(26) 



The latter equation can be rewritten for both phase-biased 
charge qubit [5J , [7] and flux qubit [[2] , taking into account 
that / = IdrcS'z, as following: I q t = I C i rc (a z ) (here a z is 
the Pauli matrix). For a phase-biased charge qubit [TT] 
the circulating current I C irc — lo is phase dependent and 
Eqs. ()24ll26p show that there are two terms contributing 
in the tank circuit's phase shift: 



tana : 



dZ 



(27) 



where Z = (a z ) is the difference between the ground and 
excited state populations. In a classical system or in the 
ground state the difference between the energy level's pop- 
ulations is constant, Z = const, and the second term in 
Eq. (|27p is zero. In contrast, for the quantum system the 
interplay of the two terms is essential, which was studied 
in Ref. [TT]. At this point it is worthwhile to notice that 
the second term can dominate at resonant excitation, as 
it was the case in the work [TTJ (c/. Figs. 3 and 5 in [TT]'). 
Hence in some cases this may be the advantage of the 
impedance measurement technique. Another advantage of 
the technique may be the possibility of the non-destructive 
measurement (see in Refs. [Tl], [5], [IS]). 

Consider now the case of a flux qubit in detail. The 
current operator is defined in the flux basis [12j , / = Ipa z , 
where Ip stands for the amplitude value of the persistent 
current, and hence the value (a z ) defines the difference 
between the probabilities of the clockwise and counter- 
clockwise current directions in the loop: (a z ) — P± — P| = 
2P l - 1. Then with Eqs. (|24II26|) we obtain 



tana « k 2 Q 



$0 0/dc' 



(28) 



where f dc = <P dc /<P - 1/2. 

For calculation of the density matrix p the Bloch equa- 
tion is conveniently used (see e.g. Ref. CIS]). This equation 
includes the relaxation and correspondingly is written in 
the energy basis. Thus, we rewrite Eq. (|28|) after intro- 
ducing the density matrix in the energy representation in 
terms of the unity matrix 1 and the Pauli matrices r^: 
p= (1/2) (l + X? x + Y? y + Z? z ) (i.e. Z is again the dif- 
ference between the ground and excited state populations) 
and obtain: 



tan a « k 2 Q 



Lip d 



( 



2A 



#o df dc \AE 



-X 



2Ip1 a f, 



AE 



Z 



(29) 



Here AE = 2y/ A 2 + (Ip^o/dc) 2 is the distance between 
the stationary energy levels and A is the tunneling am- 
plitude; about the details of this transition from current 
representation to energy representation see in Ref. [20] . 
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For one flux qubit in the ground state (X = 0, Z 
it results in the following: 



tanev 



-k 2 QL 



A 2 I 2 P 
(AE/2) S 



(30) 



It is important to notice that we obtained the result for 
the ground state, Eq. (|30"|) . which coincides with the earlier 
obtained results (see Eqs. (3-4) in [10]), but in different 
way - by differentiating the probability Pj_, Eqs. (|28H29p . 
For the description of the flux qubit in the thermal equi- 
librium one has to put X — and Z = tanh (AE/2T) in 
Eq. (|2"9")k then by plotting the phase shift versus the mag- 
netic flux, fd c , for different temperatures, one can obtain 
the suppression and widening of the zero-bias dip (that is 
in the vicinity of fdc — 0) as it was observed in the ex- 
periment presented in Ref. |21j in Fig. 3 (a), which is one 
more confirmation of our consideration. For example, the 
zero- bias dip (that is tan a at fdc = 0) is described by the 
r.h.s. of Eq. (J30j) multiplied by the factor tanh (A/T). 

If the first term in the bracket in Eq. I|29p can be ne- 
glected (which in concrete case should be checked, but this 
is usually valid for small driving amplitude <£ ac , see e.g. 
in 22J), then the expression is simplified: 



tana ~ —k 



.LI, 



8 {2I P $ f dc 



+ 



<? 
2Ip$ f dc 
AE 



dfdc 
dZ 



AE 



Of, 



dc 



(31) 



Note that at fdc — 0: a ~ Z, which means that a probes 
the changes of the upper level population. 

If a qubit is resonantly excited with the driving fre- 
quency ui, then the partial energy levels occupation prob- 
ability Z has the Lorentzian-shape dependence on fdc- It 
follows that the derivative dZjdfdc takes the shape of a 
hyperbolic- like structure, i.e. it changes from a peak to a 
dip in the point of the resonance at AE(fd c ) = hu- 



5 Conclusion 

The impedance measurement technique for the tank cir- 
cuit being coupled to the system of qubits was studied. 
The tank circuit was considered to be driven by the rf cur- 
rent and the voltage Vr to be measured. The main results 
of the work concern the phase shift a and the amplitude 
Va of the voltage Vr- It was obtained how the phase shift 
a is related to the parametric inductances of the qubits 
£i, Eq. (fl8|) . It was shown that the dynamics of the qubits 
can be studied via the amplitude as well as via the phase 
shift, Eq. (|2"2"|) . The derivations of these equations were 
presented in detail in order, first, to make all the assump- 
tions clear (small loops' inductances Lj, weak driving of 
the tank circuit Iuas , high quality factor Q and small cou- 
plings ki, slow dynamics of the tank circuit in comparison 
with qubits, u>t <C io, AE), and, second, to show how the 
parametric inductances of the qubits should be defined, 
by introducing the difference between the energy levels 



occupation probabilities, Z. We obtained that the expres- 
sion for the phase shift a in general contains both terms 
proportional to Z and proportional to dZ/dfdc, Eqs. ([27)) 
and (f3"Tjl . If the latter term dominates the resonant exci- 
tations are visualized as hyperbolic-like structures on the 
dependence of the phase shift a on the dc flux fdc- 
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